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Abstract. We discuss the effects of a weak random external potential on the 
properties of the dilute Bose gas at zero temperature. The results recently obtained 
by Huang and Meng for the depletion of the condensate and of the superfluid density 
are recovered. Results for the shift of the velocity of sound as well as for its damping 
due to collisions with the external field are presented. The damping of phonons is 
calculated also for dense superfluids. 
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1. INTRODUCTION 



The role of disorder in Bose systems has been the object of several theoretical 
works in the last years 1-7 . In this context an important issue is to understand 
the effects of disorder on the behaviour of the dilute Bose gas (DBG) where one 
expects to obtain analytic results at least for weak disorder. This problem was 
recently considered by the authors of Ref. [3] who investigated the consequences of 
a random external potential on the depletion of the condensate and of the superfluid 
component of the system. A striking result emerging from this analysis is that 
disorder is more active in reducing superfluidity than in depleting the condensate. 

The main purpose of the present work is to investigate the consequences of 
disorder on the propagation of phonons. In particular we study the shift of the 
sound velocity as well as its damping generated by scattering with the random 
external potential. 

Our formalism is based on dispersive quantum hydrodynamics at zero tem- 
perature. This formalism not only accounts in a natural way for the behaviour of 
the dilute Bose gas, but can be employed, in the hydrodynamic regime, to predict 
important properties of dense superfluids such as liquid 4 He. 

The paper is organized as follows: in Section 2 we develop the formalism of 
quantum hydrodynamics in the presence of random external potentials. In Sect. 3 
we recover, in the limit of weak disorder, the results of ref. [3] for the depletion of 
the condensate and of the superfluid density. In Sect. 4 we explore the consequences 
of disorder on the dynamic structure factor. In particular we calculate the change 
in the velocity of sound and the damping of phonons induced by disorder. The 
expression for the damping is also extended to the case of strongly interacting 
superfluids. 
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2. FORMALISM OF QUANTUM HYDRODYNAMICS. 



Let us consider a uniform superfluid system at T = 0. Our approach to the 
problem starts from a macroscopic description: we write the energy functional of 
the system in terms of the superfluid velocity, v s , and the variation of the density 
with respect to the equilibrium value, p' = p — p$. We find 8 

E=^Jd 3 r v 2 s (r) + ± J d 3 r J dV^flr - r'|)p'(r)p' (r') 

+ 1 - J d\p'(r)v 2 s (r) (1) 

The superfluid behaviour of the system fixes the irrotational law for the velocity 
field: v s = Vip. In terms of the Fourier components of the density and of the velocity 
potential, the energy functional can be written as the sum of a quadratic and a cubic 
term in the fluctuations ip^ and p' k : E = H + Hi, where 

Ho = P fJ2k 2 M 2 + 1 -J2x- 1 (k)\p^ (2) 

k k 

and 

Hi = -4= V k • k' </? k /V-k^-k' (3) 
2 ^ V k,k' 

The quantity is the Fourier transform of x -1 (r). Terms of third and higher 

order in p' have not been included in eqs. (1), (3) since in the DBG limit they do 
not contribute to the leading order. Applicability of eq. (1) is limited to the region 
of low temperatures where the whole system is superfluid. 

The quadratic functional of equation (2) gives the expression for the energy in 
the framework of the "harmonic dispersive hydrodynamics" , where the compressibil- 
ity of the system, given by is momentum dependent. In the long wavelength 
limit the hydrodynamic hamiltonian (2) coincides with the usual phonon hamilto- 
nian with the velocity of sound given by c = (po/xi®)) 1 ^ 2 - 
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Let us neglect for the moment the anharmonic term Hi and consider only 
the harmonic hamiltonian Hq. By quantizing the fields <p and p' we can cast the 
hamiltonian into the diagonal form 

H Q = y^a;fcCk g k (4) 

k 

where c k (6k) are the creation (annihilation) operators relative to an elementary 
excitation carrying impulse k and we have ignored the constant term due to the 
zero point motion. Equation (4) shows that the system can be described as a gas 
of non-interacting elementary excitations with energies 



The quasiparticle operators c k , c k are related to the operators <^ k and p k 
through the following expressions 

5^k = <£k - (0k) = (^^2) 1/2 + cl k ) 

5/5 k = /5 k -(p k )=z(|^) 1 / 2 (c k -ci k ) (6) 

where we have introduced the fluctuation operators and 5/3 k which, in the 
absence of disorder, coincide with the operators </? k and p k . 

In the dilute Bose gas limit the frequencies u>k are given by the most famous 
Bogoliubov spectrum 

u k = ^(k 2 +4p vo) 1/2 (7) 
2m 

and the compressibility x(k) takes the form 



In terms of the scattering length a the interparticle potential vq is written as: 
vo = Am a lm. 

Let us now add a perturbative field in the form of a random external potential 
U (r) coupled to the density through the interaction hamiltonian 

In eq. (9) is the Fourier transform of the external potential. In order to investi- 
gate in a simple way how the random potential affects the behaviour of the system, 
we will often make the white noise assumption in which the external potential is 
described by a single parameter R [3]: 



and where the average is assumed to be of Gaussian type. In eq. (10) m is the mass 
of the particles and the dimension of R$ is consequently (energy) 2 x (length) 3 . 

In the presence of weak disorder the expectation value of the density operator 
on the ground state can be evaluated by means of linear response theory: 

<A) = -^u 1tX (k) (ii) 

whereas the expectation value of the velocity potential operator is not affected by 
the external static field: (<£ k ) = 0. By writing both the hamiltonians H and V in 
terms of the operators Ck, cj, we get the result 

^ = -?£^i 2 >^+£^Ja (12) 

k k k 

for the harmonic hamiltonian H = H +V. The first term in (12) gives the correction 
to the ground state energy due to the external field. This quantity can be calculated 
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also beyond the dilute Bose gas approximation. In fact by applying linear reponse 
theory (eq.(ll)) to the external field (9) one finds the general result 

^ = -^E^(l^| 2 )xW (13) 
k 

holding also for strongly interacting systems. Equation (13) reduces to the first 
term in the right hand side of eq.(12) in the Bogoliubov limit where eq.(5) holds. 
It is important to notice that if one makes the white noise assumption (10), the 
correction to the ground state energy (13) contains an ultraviolet divergence since 
x(k) behaves as k~ 2 for large k. A more physical choice for y(|£/k| 2 ) would not 
yield such a divergence. 

From result (12) it turns out that, without including interaction terms among 
the quasiparticles, described by the anharmonic term Hi, the energies of the ele- 
mentary excitations are not affected by the random external field; in particular the 
velocity of sound is still given by the Bogoliubov value c = (pqVq) 1 / 2 /m. In Section 
4 we will calculate the first correction to the velocity of sound as well as its damping 
due to the random external potential. 

3. SUPERFLUID DENSITY AND DEPLETION OF THE CONDEN- 
SATE 

The current density associated with the one-fluid hamiltonian (1) is given by 

j(r) = |^=P0V s (r)+p'(r)v s (r) (14) 

where we have separeted the first and the second order terms. By taking the Fourier 
transform of the current operator we obtain 

j q = «q/w q + ^X/ k( ^ k/3 q-k ( 15 ) 
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The first term is longitudinal, whereas the second one contains both a longitudinal 
as well as a transverse component. By inserting in eq. (15) the expressions for 
the velocity potential and density (6) in terms of the quasiparticle operators c^, Ck, 
one can distinguish in a natural way between different contributions to the current 
operator. The first term in (15) is coupled to single elementary excitations: we 
call this term the one-phonon contribution to the current. The second term can be 
either the product of one quasiparticle operator and of the external field (hereafter 
referred to as disorder contribution) or the product of two quasiparticle operators 
(two-phonon contribution). 

Let us choose q in the z direction: q = (0, 0, q) and let us consider the tran- 
verse component of the current operator j^. In terms of the elementary excitation 
operators it is written as 

£ = = -1 £ i k x ) 1/2 k l 2 ( a k+a t_ h) + ff.) (16) 

where the operator jq- 2 ^ arises from two-phonon effects and does not depend upon 
the external field. 

There is an important point to stress here. The transverse current operator 
(arising from the second term in the r.h.s. of eq. (15)) is fixed by the anharmonic 
term in the hamiltonian. This implies that, for a calculation of the corresponding 
matrix elements to the lowest order, we can safely use the eigenstates and eigenval- 
ues of the harmonic hamiltonian Hq. The situation is different in the longitudinal 
channel, since the longitudinal component of the current operator (15) contains a 
term, the one-phonon contribution, which is zeroth order in the anharmonic per- 
turbation. In the next Section we will see that the anharmonic term H\ has to be 
explicitly included in the hamiltonian in order to study the longitudinal channel in 
the proper way. 
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The normal component of the fluid is obtained by considering the static trans- 
verse current-current response which, at T = 0, can be written as 

T i \ ^ I(0|J»| 2 , l(0|j-c»| 2 m . 

The normal density is then obtained from the long wavelength limit of equation 
(17) 9 

p n = \im X T (q) (18) 

q^O 

while the superfluid density is defined as p s = po — p n . In eq. (17) |0) is the ground 
state and the sum is performed over a complete set of excited states with energies 
U! n . If we now insert the expression (16) for the current operator into eq. (17) we 
can distinguish the contribution to the response function arising from two-phonon 
effects from the one fixed by disorder. The former can be shown to vanish in the 
long wavelength limit as q 2 . The latter is easily obtained by a direct calculation 
of the matrix elements using the diagonal hamiltonian (4) . The final result for the 
normal density is given by 

V k k 

The fact that the random external potential gives rise to a normal component 
in the fluid can be also regarded as the consequence of the fact that translational in- 
variance is broken by the external field and therefore the current is not conserved 10 . 

Result (19) needs some comments. First of all eq. (19) shows that the nor- 
mal component induced by the random external potential is not fixed just by the 
long wavelength properties of the system. In order to assure the convergence of 
the integral in (19) the behaviour of the elementary excitation spectrum at high 
momenta is important. A second point concerns the validity of the one-fluid hamil- 
tonian (1) when a normal component is present in the fluid. Result (19) can be 
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interpreted as the first term of the perturbative expansion of the normal density in 
the external field U. It is finally worth noticing that in the case of an ideal Bose gas 
(iVk = k 2 /2m) the integral in eq. (19) diverges at small fc's revealing the instability 
of the ideal Bose gas in the presence of disorder. 

In the DBG limit the integral in (19) can be calculated analytically by inserting 
the Bogoliubov value (7) for the energies u>k and using the white noise approximation 
for ^(|t/ k | 2 ). One gets 

P " = (a) (20) 

where n is the particle density. The applicability of result (20) is subject to the 

condition of weak disorder p n << p or, equivalently, 

m2R °«l . (21) 



Jna 

In 2D one finds 

m 3 1 

"» = g^«o- • (22) 
Both eqs.(20) and (22) are in agreement with the results of Refs. [3,11]. 

We now calculate the depletion of the condensate due to the random potential. 
To this purpose let us write the macroscopic component of the particle field operator 
\l/(r), in terms of density and phase fluctuations. The macroscopic component of 
the particle operator is formally given by 

*M(r) = V^oe^ r) (23) 

where no is the operator of the number of particles in the condensate and $(r) is the 
phase operator related to the velocity potential by the usual relation </)(r) = $(r)/m. 
In the DBG limit, where no — n = po/m, the following expansion is valid up to 
first order in the density and phase fluctuations 

* M (r) = v^(l+2^/(r) + i£(r)) ( 24 ) 
10 



By taking the Fourier transform of eq. (24) one finds 

5a q = a q - (a q ) = iyfiiSQ^ + — U^p q (25) 

with (d q ) = v^q.o + l/(2v / ^- m )(P q ) (we have used ($ q ) = 0). The momentum 
distribution is thus given by 

n{q) = (Sa^aJ + |(a q )| 2 = ^<<5*- q **Q> + " \ + 

(26) 

The first term in the right hand side of equation (26) represents the contri- 
bution to the momentum distribution given by the phase fluctuations. In the long 
wavelength limit this term gives the main contribution and is responsible for the 
well known 1/q divergence of n(q) 12 . In the region of higher momenta also the 
density fluctuations (5/9_ q 5/3 q ) become important. The constant term —1/2 arises 
from the commutation relation between density and phase operators. It is worth 
noticing that the random potential explicitly affects the result for the momentum 
distribution because the expectation value of the particle operator d q on the ground 
state is no longer zero for q ^ 0: 

«) = 7w^<P q > = -Tr^^ff^xiq) (27) 

By using the Bogoliubov results for the density and phase fluctuations (correc- 
tions to these quantities due to disorder give rise to higher order effects) we finally 
obtain the result 



^) = ^(^ + -o-^) + ^ q | 2 >g (28) 



4 

where the first term coincides with the prediction of the Bogoliubov model. The 
depletion is obtained by integrating result (28) and we find the usual Bogoliubov 
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value with a correction proportional to p n : 

An = n-n = ^ = ^ {naf 2 + ^p n (29) 

Result (29) is in agreement with the findings of Ref. [3]. In particular it shows 
that the relative depletion of the condensate due to disorder is a factor 3/4 smaller 
than the corresponding reduction of the superfluid density. This result holds for 
any choice of the average y(\Uk\ 2 ) since the same integral over momenta is involved 
in the calculation of both p n and An (see eqs.(19) and (28)). We stress however 
that eq. (29) holds only in the presence of weak disorder. It cannot be easily 
extrapolated to the large external fields. 

4. VELOCITY AND DAMPING OF SOUND 

In the present Section we discuss the effect of disorder on the dispersion of the 
phonon mode at T = 0. To this purpose we first calculate the compressibility of 
the system, x(0), defined by the relation 

r) 2 F 

x- 1 ® = Wv (30) 

where E is the ground state energy of the system. The contribution to E due to 
disorder is given by eq.(13) or equivalently, in the Bogoliubov limit, by the first 
term in the r.h.s. of eq. (12). After taking the second derivative with respect to 
the density we find the result 

1 k 8 

X-\0) = Xo\0)(l + y d^l 2 )^) (31) 

k fc 

where Xo(0) is the compressibility of the DBG. By making the white noise assump- 
tion (10) it is possible to relate result (31) to the normal density p n given by eq.(19). 
We find 

X - 1 (0) = Xo 1 (0)(l + ^) (32) 
12 



It is worth noticing that, differently from the correction to the total energy as 
well as to the first derivative dE/dp (proportional to the chemical potential), the 
contribution of disorder to the compressibility is well defined also within the white 
noise assumption. This is due to the fact that the quantity d 2 x(k)/dp 2 , differently 
from x and dx(k)/dp decreases as k~ A when k — > oo and the corresponding integral 
is well behaved for large k. 

We are now ready to calculate the velocity of sound. Using the relation 

c 2 = -t^t ( 33 ) 

x(o) 1 ) 

typical of superfluids and eqs. (20) and (32) for p n and x(0) we obtain the result 

^ = < 1 + i!) (34) 

holding within the white noise assumption. Equation (34) shows that the velocity 
of sound increases with disorder and that consequently the phonon peak in the 
dynamic structure function 

/+oo 
dte^([5p_ q (t),5p q (0)]) (35) 
-oo 

is shifted to the right according to the law 

S phonon (q,uj) = ^5(uj - cq) . (36) 

The normalization factor p s q/2c ensures that the phonon exhausts the compress- 
ibility sum rule 

X (0) = 2lim q ^ [ -S(q,oj)dw (37) 

consistently with relation (33) for the sound velocity. 
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In the second part of the section we calculate the damping of phonons due 
to collisions with the external potential. To this purpose we calculate directly the 
propagator defined by the time ordered product 

Dn^t) = -z(T(p q (t)p'_ q (0))> = -z((0(t)p q (t)p'_ q (O) + 0H)p'_ q (O)p q (t))> (38) 

in terms of which the dynamic structure function is written as S(q, u) = —^ImDn(q, 
As anticipated in Section 2 in order to calculate to the proper order the effects 

of the external random potential in the longitudinal channel it is essential to include 

the anharmonic term (3) in the hamiltonian of the system. 

We calculate the imaginary part of the propagator Dn(q,u) perturbatively in 

the interaction hamiltonian 

k,k' k 

up to terms of second order in the random potential V. The calculation is straight- 
forward and the details are given in the Appendix. The relevant contribution to the 
imaginary part of the propagator arises from the diagram of Figure 1. By adding 
this contribution to the phonon peak (36) we find, for small q, the result 

S{ q , U ) = %Hu-ct) + ^ ^ g2)2 § £ ±{W)H%sS{u - Mt ) (40) 

where the second term holds for u> not too close to cq: \u — cq\ > cq. In the 
hydrodynamic limit (q — > 0, u — > 0) result (40) becomes 

Q( \ P&Xt \ Q 20j5 P° 1 l\TT |2\ (A-i\ 

S(q,u)=—S(u-cq) + {oj2 _ c2q2)2 127r 2 c7 y<l^l >H4=«-/c (41) 

Eq. (41) corresponds, in the limit of weak disorder (see eq.(21)), to the first two 
terms of the series of the most divergent diagrams characterizing the behaviour of 
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the propagator Dn, near the pole. The sum of this series yields the following result 
for S(q,u>) near the pole 

q( ) = m! 2c Q r (g) U2) 

[q ' UJ> tv (u 2 -c 2 q 2 ) 2 + 4c 2 q 2 T 2 (q) [ ' 

with the width T(q) given by 

r(9) = i?<4 (|£7 - |2> ■ (43) 

Eqs. (42), (43) explicitly show that the phonon peak is broadened by the presence 
of disorder. 

It is not difficult to extend results (42), (43) to the case of a strongly interacting 
Bose superfluid where quantum hydrodynamics is expected to provide the proper 
description in the hydrodynamic regime. To this purpose we have to add to the 
anharmonic term (3) a contribution cubic in the density fluctuations yielding, for 
Hi, the more general expression 13 

1 1 / d c 2 \ 

Hi = — f= ^ k ' k'^kPk'-k^-k' H f= I i lY^ PkPk'-kP-k' ( 44 ) 

aWKdppJ^ 

By accounting for the diagrams relevant for the calculation of the damping we get 
the relevant result (see the Appendix for details) 

r W = i?F«^( 1 + 3 ?^7 )2 ) ' (46) 

In superfluid 4 He the second term in eq. (45) turns out to be of the same order as the 
first one. Vice-versa in the DBG the derivative with respect to the density of c 2 / p 
is of higher order in the scattering length and its contribution can be consequently 
neglected. Notice that the damping mechanism discussed above is dominant at 
small q compared to the damping r p h_ p h due to anharmonic interactions among 
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phonons, which exhibits a q 5 law (see Ref. [13]). In the dilute Bose gas and using 
the white noise approximation we find the result 

m = ™^i_ (46) 

T ph - ph (q) 3 po qa 

Starting from result (40) we can finally calculate the low q behaviour of various 
energy moments of the dynamic structure factor. These moments can be written 
as the sum of the phonon contribution (first term in eq. (40)) and of the collisional 
term (second term of eq. (40)). The latter contribution is always proportional to 
q 2 . For example, by integrating eq.(40) with respect to u>, one finds the following 
result for the low q expansion of the non energy weigthed sum rule (static structure 
factor): 

Sta) = jf*-Sta, W )=M + 9 ^/ A _I_^)M|! (47 ) 

By making the white noise approximation for y(\Uk\ 2 ), one can relate the second 
term to the normal density (19): 

s{q) = ^ + q Pn m^Ta (48) 

Eq. (48) shows that the static structure factor, at low q, is exhausted by the phonon 
peak. One immediately shows that the phonon peak exhausts also the inverse energy 
weighted sum- rule fixed by the compressibility of the system (eq. (37)). 
For the energy weighted sum rule we find the result 



f°° q 2 q 2 q 

m\ = J duu)S{q,u) = p s — + p n — = p — (49) 



which coincides, as expected, with the model independent f-sum rule, proportional 
to the total density of the system. In this case the contribution of the collisional 
term is crucial in order to satisfy the sum rule also at low q. This contribution 
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can be always expressed in terms of the normal density, independently of the white 
noise assumption. It is interesting to remark that a similar behaviour is exhibited 
by transverse spin excitations (magnons) in antiferromagnets 14 . Also in this case 
the magnon exhausts only a fraction of the energy weighted sum rule, the remaining 
part being, in this case, exhausted by multimagnon excitations. The one-magnon 
contribution to the energy weighted sum rule is proportional to the spin stiffness 
coefficient which plays the role of the superfluid density. This analogy between 
disordered bosons and antiferromagnets is due to the fact that in both cases the 
current is not conserved due to the lack of translational invariance. 

The fact that the phonon peak does not exhaust the f-sum rule implies that 
the Feynman approximation 

= m> (B0) 

for the energy of elementary excitations does not coincide with the phonon disper- 
sion at low q. Actually using result (47) for S(q) one finds ujf(q) — > cqpo/p s a value 
higher than cq. It is also worth noticing that the Feynman ratio (50) is affected by 
the q 2 correction in S(q). Such a correction is absent in translationally invariant 
systems. 
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APPENDIX 



In this appendix we evaluate the imaginary part of the density-density prop- 
agator Dn(q, uj) in the long wavelength limit, to first order in R . We use time 
dependent perturbation theory with the perturbative field given in equation (39). 

We need the free field propagators: 

/+oo 2 
dt e^(T(^(t)p'_ q (0))) = i 2 ™ 2 ,. n (Al) 

/+oo -i 
dt e™*<r(0 q (f)0_ q (O))) o = t-SL-5 ^— - (A2) 
-oo POQ ^ - + «u 

/ + oo 
eft e^(T(p q (t)^_ q (0))) o = - 2 _o , - n (A3) 
-OO W ([ tU 

where the subscript after the brackets (...)o means that the expectation value is 
taken on the ground state relative to the unperturbed hamiltonian Hq. 

The density-density propagator is obtained from the general formula of pertur- 
bation theory 

iDni^t) = ^<T(p q (t)p'_ q (0)S)> (A4) 
where the time evolution operator S is defined by the series expansion 

00 ( —i) n f+ oc r+oo 

S = J2—T dt!... dt n T(H int (t 1 )...H int (t n )) (A5) 

n=0 n - J ~°° J -°° 

It is better to use diagrams to represent the different terms arising from equa- 
tion (A4) after using Wick's theorem for time ordered products. The anharmonic 
potential Hi is represented by a vertex with three lines, whereas to the external 
random potential corresponds a vertex with just one line. 

At finite frequencies the leading correction to the imaginary part of Dn(q,w) 
comes from terms which are fourth order in the interaction hamiltonian Hi nt . These 
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terms correspond to diagrams containing two vertices relative to the anharmonic 
potential and two random potential vertices. 

It turns out that only one diagram is relevant and this is given in figure 1, 
where we have represented by a full line, D 22 by a dashed line and D® 2 by a 
long-dashed line. In fact the external random field is time independent. This means 
that only the density-density propagator, which has a zero frequency component, 
can be connected to the random potential vertices. The contribution arising from 
the diagram in figure 1 to the density-density propagator is given by 

D 11 (q 1 oj) = D° 11 (q,oj) + -^(D^uj)) 2 

k 

The square of the free propagator in equation (A6) can be approximated, for 
u> not too close to the pole (\u — u q \ > u q ), with the expression 

(OU^i^jJ^ (A7) 

Result (40) for the dynamic structure factor follows directly from equations (A6) 
by using the relation 5'(q, u) = —-^ImDn(q,u) and after taking the q — > limit. 

The imaginary part of the density-density self energy can be directly evaluated 
from equation (A6). In the long wavelength limit one gets 

-JmSute, uj) = I -^1 U\U^\ 2 )kl^5{u - u k ) (A8) 
z PoQ v k k 

Result (43) for the damping in the phonon region follows from eq. (A8) by using 
the relation 

r( ? ) = M(-/mEn(g,a; = c ? )) (A9) 
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In the case of strongly correlated systems the damping of the phonon mode 
can be obtained in the long wavelength limit by using the anharmonic term (44). 
In this case another diagram is relevant to the calculation of the damping and it is 
given in figure 2. The imaginary part of the self energy is given by 

_ /mSll(9 , u) = £1^ i(l^) + <f ) 2 5) (Aio) 

and result (45) follows directly from relation (A9). 

The present diagrammatic technique could be used for the calculation of higher 
order corrections in the parameter Rq. To perform the corresponding averaging one 
must assume gaussian statistics for the random external field. 
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FIGURE CAPTIONS 



Fig.l. The relevant diagram of the fourth order in H int (eq. (39)). The full line 
represents D® 1: the dashed line D 22 and the long-dashed line D® 2 . 
Fig. 2. The other relevant diagram of the fourth order in Hint = Hi + V with Hi 
given by eq. (44). The full line represents D^. 
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